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Abstract. Studying the internal structure of extragalactic jets is crucial for understanding their physics. The 
Japanese-led space VLBI project VSOP has presented an opportunity for such studies, by reaching baseline 
lengths of up to 36,000 km and resolving structures down to an angular size of as 0.3 mas at 5 GHz. VSOP 
observations of the jet in 0836+710 at 1.6 and 5 GHz have enabled tracing of the radial structure of the flow on 
scales from 2 mas to 200 mas along the jet and determination of the wavelengths of individual oscillatory modes 
responsible for the formation of the structure observed. We apply linear stability analysis to identify the oscillatory 
modes with modes of Kelvin-Helmholtz instability that match the wavelengths of the structures observed. We 
find that the jet structure in 0836+710 can be reproduced by the helical surface mode and a combination of the 
helical and elliptic body modes of Kelvin-Helmholtz instability. Our results indicate that the jet is substantially 
stratified and different modes of the instability grow inside the jet at different distances to the jet axis. The helical 
surface mode can be driven externally, and we discuss the implications of the driving frequency on the physics of 
the active nucleus in 0836+710. 
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1. Introduction 

Resolving the transversal (or radial) structure in extra- 
galactic jets is a crucial step in our understanding of the 
physics of these objects. On parsec scales, this has become 
feasible only recently, using spac e observations 
with VSOF0 (jLobanov fc Zensus I [iuOlh . These observa- 
tions revealed the presence of a double helical structure 
inside the jet of 3C 273, which can be attributed to a com- 
bination of helical and elliptic modes of Kelvin-Helmholtz 
(KH) instability. N umerical simulations further support 
this interpretation ( Perucho et al. 20061 ). In this letter, 



we expand the scope of this investigation and use VSOP 
observations of the radio jet in the quasar S5 0836+710 at 
1.6 and 5 GHz to estimate basic physical properties of the 
relativistic flow, focusing specifically on the radial profiles 
of its velocity and density. 

The luminous quasar S5 0836+710 at a redshift z = 
2.16 hosts a powerful radio jet extending up to kilopar- 
sec scal es (iHummel et al. 1992 ). VLBI monitoring of the 



source (jOtterbein et al. 1 119981 ) has enabled estimates of 
the bulk Lorentz factor 7j = 12 and the viewing angle 
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9j = 3° of the flow. Presence of instability developing in 
the jet is suggested by the kink structur es observed on mil 



liarcs econd scales with ground VLBI (jKrichbaum et al. 
19901) . 

In the VSOP image of 0836+71 at 5 GHz, oscillation s 



of the jet ridge line are observed (jLobanov et al. 1 119981 ). 



with wavelengths of 7.7 mas and 4.6 mas. These structures 
have been identified with the helical surface and ellip- 
tic surface modes of KH instability, yielding estimates of 
basic physical properties of the flow: the Lorentz factor 



7j 



11, Mach number Mj ~ 6 and jet/ambient medium 



density ratio f?(pj/p a ) = 0.04. High dynamic range VLB^O 
observations of 0836+710 at 1.6 GHz show the presence 
of an oscillation at a wavelength as long as ~ 100 mas 
( Lobanov et al. I l2006l) . This wavelength cannot be recon- 
ciled with the jet parameters determined from the two 
shorter wavelength oscillations, indicating that the flow 
may have a complex, stratified transversal structure in 
which emission at lower frequencies originates from outer 
layers of the flow. 

Taking into accou nt that the parameters given in 
Lobanov et al. I (jl998l ) were obtained under the assump- 
tion of a vortex sheet transversal structure of the jet, we 
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explore the possibility that the jet is sheared, with the 
long wavelength corresponding to a surface mode growing 
in the outer layers, whereas the shorter wavelengths are 
identified with body modes developing in the inner radii of 
the jet. In our approach, we assume that the magnetic field 
is not dynamically important and that the structure of the 
jet is due to KH instability. Under these assumptions, the 
jet structure and evolution can be described within the 
framework of linear stability theory for cylindric outflows. 

In Section 2, we describe the method used to solve the 
linear stability problem for cylindric relativistic flows and 
provide the respective so l utions for the set of parameters 
given in lLobanov et al7 (1998). This solutions are com- 
pared in Section 3 with the wavelengths observed in the 
jet of 0836+710. Section 4 summarizes the main results 
of this work and puts them in the broader context of the 
physics of relativistic outflows. 

2. Linear analysis 

We describe the flow in cylindrical coordinates (r, z, <fi), 
with r and z defining the radial and axial directions, re- 
spectively. We consider a sheared transition in the axial 
velocity, v z and the rest mass density, p, between the 
jet and the ambient medium of the following form (see 
I Perucho et al. 1120071 for a detailed description of the linear 
analysis) a(r) = + (ao — a^)/ cosh(r m ), where a(x) is 
the profiled quantity {v z or p) and ao, are its values at 
the jet symmetry plane (at r = 0) and at r — > 00, respec- 
tively. The integer m controls the steepness of the shear 
layer. In the limit m — > 00, the co nfiguration turn s into the 
vortex-sheet described bv lHardeel (|2000h . The jet 

and the ambient medium are in pressure equilibrium. An 
adiabatic perturbation of the form oc <?(r) exp[i(k z z — tot)] 
is introduced in the equations of the flow. Here, u and k z 
are the frequency and longitudinal wavenumber of the per- 
turbation, respectively, and the function g(r) defines the 
radial structure of the perturbation. The resulting pres- 
sure perturbation P\ is then obtained from a second order 



ordin ary differential equation (equation 13 in lBirkinshaw 



Il984h . This equation can be solved by integrating from 
the jet axis to a point outside the jet, where boundary 
conditions on the amplitude of pressure perturbation and 
its first derivative are imposed. We apply the following 
boundary conditions: a) symmetry or antisymmetry of the 
perturbation on the jet axis, b) the Sommerfeld condition 
requiring the solutions to approach zero at infinity, and 
c) no incoming waves at infinity. These conditions ensure 
that the solutions can be given by Bessel functions inside 
the jet and Hankel functions outside. The solutions are 
obtained in the spatial domain, i.e., assuming uj is real 
and k is complex. In this description, 1/I(k) gives the 
growth length (or e- folding length), defined as the dis- 
tance at which the perturbation amplitude increases by 
one exp onential factor. We adopt th e jet parameters de- 



rived in Lobanov et al. I (|1998L [2006) and solve the pres 



sure perturbation equation with several different values of 
m. The resulting solutions are plotted in Fig. Q] for the 



helical instability modes and for two values of to: to = 8 
(left panels) and to = 200 (right panels). The m = 8 case 
corresponds to a broad shear layer of ~ 0.6 Rj in width, 
whereas the to — 200 one implies a shear layer width of 
~ 0.1 Rj, i.e., approaching the vortex-sheet case. 

The plots in Fig. [1] indicate that, in thicker shear lay- 
ers (m — 8), the surface mode grows faster than the 
body modes. The minimum growth lengths of the surface 
mode indicated by vertical dashed lines in Fig.[T]) are real- 
ized at long wavelengths (horizontal dashed lines show the 
wavenumber of the mode): Ag H 8 ~ 160 Rj. The smallest 
growth lengths of all body modes is achieved by the fifth 
and sixth order body modes at Ab6,8 ~ 14 -Rj) whereas the 
minimum growth length of the first body mode (expected 
to be among the most prominent in the jet) is realized at 
Abi >8 ~40i?j. 

For a thin shear layer (m = 200), the solutions ap- 
proach the vortex-sheet approximation at the lowest fre- 
quencies (longest wavelengths). The minimum growth 
lengths of low order helical body modes are similar to 
that of the helical surface mode, with the minimum real- 
ized for the first body mode at Abi,200 ~ 25 -Rj . The min- 
imum growth length of the surface modes shifts to higher 
frequencies and shorter wavelengths (Ag H 200 ~ 100 -Rj) 
compared to the respective values obtained for to = 8. 
For to < 8, the growths of all the modes are strongly re- 
duced as the width of the shear layer is increased. For 
to > 200 the solutions converge to the vortex-sheet case 
also at the highest frequencies/shortest wavelengths. 

3. Results 

The solutions of the linear stability problem obtained in 
the previous section can be compared with the wave- 
lengths observed in the jet in 0836+710. This compar- 
ison requires an estimate of the radius of the jet, in 
order to convert the frequencies and wavenumbers of 
the solution into physical units. The true jet radius can 
be estimated from the VL BA radio maps presented in 
Lobanov et al. ] (|l998l . I2006T) . using the following relation 
Rj(mas) = Q.hy/Df — 6 2 , where is the observed width 
of the jet and b is the beam width transversally to the 
direction of the jet. We m easure the diameter of the jet at 
1% of the peak emission dWehrle et al. Ill992h at the base 
of the jet, implicitly assuming that the outer parts of the 
shear layer emit less than 1% of the radio power at these 
frequencies. This yields the jet radii of Rj,i.6 ~ 17mas 
and Rj t 5 ~ 0.64 mas at 1.6 GHz and 5 GHz, respectively. 
With these values, the intrinsic (rest frame) wavelengths, 
A™° d , can be obtained and expressed in units of Rj. In or- 
der to compare these wavelengths with the observations, 
the observed wavelengths must be corrected for relativis- 
tic motion, projection and cosmological effects, yielding: 
Amt = A obs (1 + z) (1 - v w /c cos0j)/ sin#j, where A obs is 
the observed wavelength, z is the redshift, v w is the wave 
speed and 0j is the jet angle to the line of sight. The wave 
speed v w — u/k is obtained from the solutions to the lin- 
ear problem. The wave speed of the helical surface mode 
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Fig. 1. Real (upper curves) and imaginary (lower curves) parts of the wavenumber versus frequency for the helical 
modes and the parameters given for the jet in 0836+710. The vertical dashed lines indicate the minimum growth 
lengths of the surface modes, whereas the vertical dotted lines indicate the minimum growth length of the first body 
modes for m — 200 and the minimum growth length of all body modes for m = 8. The horizontal lines indicate the 
wave number at which the minimum growth length occurs. S stands for the surface mode, Bl for the first body mode, 
and B5-B6 indicate the body mode that gives the minimum growth lengths for m = 8. 

Table 1. Intrinsic wavelengths corresponding to the observed structures in the jet of 0836+710 





Aint(100 mas) 


Ai n t(100 mas) 


A int (7.7 mas) 


A int (7.7 mas) 


A int (7.7 mas) 




(i)w = 0.6 c) 


(« w = 0.8 c) 


(v w = 0.2 c) 


(». = 0.5 c) 


(v w =0.7c) 


Rj, 1.6(17 mas) 


140i?j 


72 Rj 


22 Rj 


14 -Rj 


SRj 


i?j,5(0.6 mas) 


4 10 3 Rj 


2 10 3 Rj 


620 Rj 


390 Rj 


230 Rj 



Notes: The wave speeds are given at the minimum growth length of the relevant modes: the helical surface mode for m = 8 
(f w = 0.6 c) and m = 200 (w w = 0.8 c), the fifth and sixth body mode for m = 8 (w w = 0.2 c), and the first body mode for m = 8 
(«w = 0.5c) and m = 200 (« w = 0.7c). 



(most likely responsible for the longest observed wave- 
length of 100 mas) ranges from v w ~ 0.6 c to v w ~ 0.8 c, 
depending on the thickness of the shear layer. The dom- 
inant body modes (responsible for the shorter observed 
wavelengths of 7.7 mas and 4.6 mas) have speeds ranging 
from v w ~ 0.2 c to v w ~ 0.7 c. Table [T] gives the intrinsic 
wavelengths for two main observed wavelengths (100 mas 
and 7.7 mas), calculated with different appropriate wave 
speeds in units of the jet radius. The respective intrin- 
sic wavelengths of the 4.6 mas mode can be obtained by 
dividing the values for the 7.7 mas mode by a factor of 
1.7. 

We note that the wavelengths obtained for the m = 
200 case cannot be reproduced by any combination of 
excited KH modes at, or close to, their resonant fre- 
quencies, although the lack of error estimations prevents 
us from completely ruling out this possibility. This fact 
alone points towards the jet being transversally stratified. 
Further evidence for this comes from a detailed analysis 
of the broad shear layer case (to = 8). 

For the m = 8 case and for i?j,5, the intrinsic wave- 
lengths for the longest structure (100 mas) have no coun- 
terparts in the solutions to the linear problem, as already 
pointed out bv lLobanov et al.l ( 20061 ). This suggests that 



the radiating particles emitting at 5 GHz occupy a cen- 
tral spine in a stratified jet, while the longer-wavelength 
structure seen at 1.6 GHz is generated in outer layers of 
the flow. 

The intrinsic wavelengths obtained for i?j,i.6 are close 
to the most unstable wavelengths provided by the to = 
8 shear layer. In this case, the helical 100 mas (140 Rj) 
structure, corresponding to a wavenumber k ~ 0.045, can 
be identified with A* H - 160 Rj (k* K ~ 0.04). The helical 
7.7 mas (22 structure can be identified with the fifth 



or sixth body mode excited at A > Ag 6 H ~ 14i?j, and 
the elliptic 4.6 mas (13i?j) structure would correspond to 
the fifth and sixth elliptic body modes growing at Ag 6 E ~ 
14 Rj, as indicated by the linear solution (not shown). We 
consider this a better explanation, and we adopt -Rj,i.6 as 
the generic jet radius Rj in the following discussion. 

One possible alternative to the stratified jet scenario 
is to consider de viations of the b a sic jet param eters from 
those derived in 
this alternative, 



Lobanov et al~l (|l998l l2006h . To assess 



we have calculated the solutions of the 
stability equation for different sets of parameters, varying 
the jet Lorentz factor, rest-mass density ratio and spe- 
cific internal energy. Our conclusions are summarized be- 
low: 1 ) We find that for colder jets (with a sound speed 
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c s j ~ 0.054, corresponding to M ~ 18) the growth lengths 
(k^ 1 > 10 3 Rj) of the surface and first body modes are 
too large to explain the growth of the instability in the 
jet. 2) For hotter jets (c s j ~ 0.55, M ~ 1.8), no body 
modes exist that could reproduce the observed 7.7 mas 
and 4.6 mas wavelengths. At the same time, the surface 
modes alone cannot reproduce the three wavelengths ob- 
served. 3) Lighter jets (rj = 4 ■ 10~ 4 ) produce solutions 
similar to those arising for colder jets, with the excep- 
tion that the growth lengths are slightly shorter in this 
case, but still of the order of 500 Rj, i.e., even longer than 
the longest observed wavelength taking the radius of the 
jet at 1.6 GHz. 4) Denser jets (77 = 0.4) show the sur- 
face and first body modes at the appropriate wavelengths, 
but the surface mode grows at a too long growth length 
(fc.r 1 ~ 500 Rj). Moreover, the observations indicate that 
the surface mode must grow faster than the body modes, 
contrary to the results derived for a dense jet. 5) Slower 
jets (7j = 2) show very short growth lengths for the sur- 
face mode (fc" 1 ~ 17Rj) and the solutions yield too short 
wavelengths for the body modes. 6) Finally, much faster 
jets (7j = 100) produce too long growth lengths. This 
analysis indicates that the jet in 0836+710 must have a 
Lorentz factor close to 12 or slightly smaller, a density ra- 
tio of 10 _2 -10 _1 and a sound speed c s j ~ 0.2 or slightly 
la rger — not deviating much fr om the parameters derived 



Lobanov et al~l (|l998l . l2006th 



4. Discussion 

The description developed in this paper explains all three 
major oscillations observed in the jet of 0836+710, if the 
flow speed and density are both transversely stratified. 
The longest observed wavelength of 100 mas represents 
the helical surface mode, whereas the two shorter ones 
(7.7 mas and 4.6 mas) are identified with the fifth/sixth 
order helical and elliptic body modes, respectively. It is 
also possible that the short wavelengths correspond to a 
lower order body mode (first body mode, for example) 
growing not at a resona nt wavelength. 

It has been shown bv I Hardee et al. I ( 1994)that the he- 
lical surface mode of KH instability can be driven by an ex- 
ternal periodic process. In this case, coupling this process 
to the helical surface mode requires the driving frequency 
to be lower than the resonant frequency. For 0836+710, 
the driving frequency of the helical surface mode at its 
minimum growth length is 0.025 c/i?j would imply a driv- 
ing period of Td r ~ 5.6 • 10 7 yrs. This is simi lar to the driv- 
ing p eriod of ~ 2 • 10 7 yrs found for 3C 449 (jHardee et al. 
19941) . 



The period T<j r can be produced by a number of phys- 
ical processes, including a binary bla ck hole and a mis - 



aligned torus as the most plausible (lAppl et al. I 1996) 



The latter can however be discarded after iLu & Zhou 



( 20051) iave shown that the absolute magnitude of the 
galaxy in the B band anticorrelates with the precession pe- 
riods due to misalignment between the black hole rotation 
axis and the outer regions of the accretion disk. The abso- 



lute magnitude of 0836+ 710 in the R band is M R = -29.9 
( Caccianiga et al. 1 [2002) . Assuming a similar magnitude 
for the B-band yields precession periods of 10 3 — 10 5 yrs, 
depending on the a ccretion disk and black-hole properties 
(ILu & Zhoull2005h . This is much shorter than Td r we have 
derived. Combining Tdr with the es timate of the black hole 
mass in 0836+710 (2-10 8 -10 9 M ? ; lTavecchio et al. II2OO0I) 
yields a companion mass of 10 4 — 10 7 M , depending on the 
separation between the two black holes (10 17 — 10 18 cm). 

Also, the kiloparsec-scale structure of 0836+710 shows 
a large, decollimated feature ( Hummel et al. |[l992l ). with- 
out any emission between it and the 1.6 GHz jet. It re- 
mains to be tested by new high dynamic range observa- 
tions whether the long instability mode ultimately causes 
the disruption of the flow, and produces this decollima- 
tion. 
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